Experimental Realization of the Deutsch- Jozsa Algorithm with a Six-Qubit Cluster State 
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We describe the first experimental realization of the Deutsch- Jozsa quantum algorithm to evaluate the prop- 
erties of a 2-bit boolean function in the framework of one-way quantum computation. For this purpose a novel 
two-photon six-qubit cluster state was engineered. Its peculiar topological structure is the basis of the origi- 
nal measurement pattern allowing the algorithm realization. The good agreement of the experimental results 
with the theoretical predictions, obtained at ~lkHz success rate, demonstrate the correct implementation of the 
algorithm. 
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PACS numbers: 03.67.Ac 03.67.Bg 03.67.Mn 

Introduction. - In the last decade, quantum information 
processing and, in particular, quantum computation, have 
been conquering increasing interest and importance in the 
scientific community, supported by the promising theoretical 
and experimental results obtained. One of the many present 
efforts is the construction of quantum hardware, which up to 
now has been realized by following different experimental 
techniques IHI2]. In this way, it was then possible to demon- 
strate the correct functioning of one and two-qubit logic 
gates as well as the successful implementation of quantum 
algorithms which strongly show the efficiency of a quan- 
tum computer with respect to its classical analogue. Among 
these, the Deutsch- Jozsa (DJ) algorithm is the first example 
of the speed-up exhibited by a computer taking advantage of 
quantum mechanics in the evaluation of a global property of 
an n-bit boolean function ||3|]. 

In this Letter we report the realization of the Deutsch- 
Jozsa algorithm in the framework of the one-way model 
of quantum computation which has already proved 

successful in the construction ofquantum gates such as the 
controlled-NOT (CNOT) gate 041 and in the implemen- 
tation of the Grover 10, @- U ] and the Deutsch algorithms 
fUHH]. The latter corresponds to the case n = 1 and is based 
on the use of four-qubit cluster photon states. Here we get 
the access to the case n = 2 by taking advantage of a pecu- 
liar two-dimensional two-photon six-qubit cluster state gen- 
erated by a source of multi-qubit cluster states whose perfor- 
mances have been already demonstrated (8, 13]. At variance 
with the simple case n = 1, the DJ algorithm allows to take 
advantage of the exponential growing of the computational 
speed-up for increasing values of n, as said. Hence the re- 
sults presented in this paper are important in that they open 
the way to the implementation of the DJ algorithm with a 
still larger number of qubits. Although the DJ algorithm has 
been implemented before with photons 11411 . our realization 
represents the first realization with a 2-bit function in the 
context of measurement-based quantum computation. 
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FIG. 1: Graphs associtated to (a) the |E) cluster and (b) the HEb) 
hyperentangled state. Qubits 1, 4 are encoded in the E/I momen- 
tum, qubits 2, 5 are encoded in polarization (H/V) and qubits 3, 6 in 
the i/£ momentum, (c) Annular section of the conical SPDC emis- 
sion of a Type I phase matched crystal. The source produces the 
|HEe) state over eight spatial modes. Two half wave-plates (A/2) 
with vertical optical axis intercepting three modes of the A photon 
are used to transform |HEe) into the |E) state. 



Realization pattern for the Deutsch-Jozsa algorithm. - 
Let us briefly recall the generalized version of the Deutsch 
algorithm yfl, where a set of n qubits constitutes the input of 
a black box, usually known as the Oracle, which implements 
the n-bit boolean function f(x) such that / : {0, 1}™ — > 
{0, 1}. The aim of the DJ algorithm is to determine whether 
the function evaluated by the oracle is constant or balanced; 
a function / is said to be balanced if it is equal to when 
calculated in half of the possible values of x and equal to 
1 when the remaining allowed values for x are taken into 
account. Classically, 2"~ 1 + 1 queries to the oracle are nec- 
essary to solve the problem while, in the frame of quantum 
mechanics, the answer comes with one single query. Con- 
sequently, the greater is the number n of qubits involved, 
the more evident is the difference in the performances of the 
quantum computer with respect to its classical counterpart. 
The initial state of the system is |0) ® |0) ® • • • <g> |0) = |0)® n 
and an ancillary qubit in the state |1) is added to the n input 
qubits. The operation performed by the oracle is given by 
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FIG. 2: Balanced function: (a) algorithm realized according to the 
pattern of single-qubit measurements and (b) equivalent circuit im- 
plementing the DJ algorithm for n — 2. Gray gates represent Pauli 
errors. 



FIG. 3: Constant function: (a) algorithm realized according to the 
pattern of single-qubit measurements and (b) equivalent circuit im- 
plementing the DJ algorithm for n = 2. Gray gates represent Pauli 
errors. 



(H <g> H)[//(l-r™ ® H), where H is the Hadamard gate; 
the unitary operator U / acts on the states of the computa- 
tional basis so that Uj\x) \y) = \x}\y(B f{x)}. The final state 

is found to be £;^(-l) /(8)+x ' v |y>)|l}. Measuring 
the state of the n qubits in the computational basis leads to 
the conclusion: if we get the state |0)® n the function / is 
constant, otherwise it is balanced, as seen from the above 
expression for the final state of the global system. Moreover, 
the measurement of the ancillary qubit in the computational 
basis is expected to always return the |1) state. 

We now go into the details of the proposed experimen- 
tal realization of the DJ algorithm for a function acting on 
n = 2 bits. In this case, the boolean function which we are 
interested in is such that /: {0,1} 2 — s-{0,l}. The function 
/ can be calculated in its four arguments x = 0, 1, 2, 3, with 
x = 2x\ + xq and xq,x\ = 0, 1. Among the 16 possible 
functions of this kind, we focus our attention on the bal- 
anced function such that /b(0) = /s(3) = 0, /s(l) = 
/s(2) = 1, and on the constant function fc for which we 
have that fc{%) = for every allowed value of x. We thus 
identify the state \x)q = \xi)-\xq)2 as the input entering the 
oracle. In the former expression the subscripts 1 and 2 refer 
to logical qubits. As we know, the implementation of the DJ 
algorithm requires an additional ancillary qubit in the initial 
state \y) a = \y)3> where 3 is the logical qubit associated 
to the ancilla. For the previously defined functions we have 



that Uf c = 11 and U/ B = CNOTjgCNOT^, respectively 1 . 

In the framework of one-way quantum computing, the 
starting point of any computation is the construction of a 
multi-qubit cluster state; successively, the choice of a se- 
quence of single-qubit measurements determines the pro- 
gram to be executed on the quantum computer. For a re- 
view on graph and cluster states and their use for one-way 
computing see J3, El H El ■ Let us start from the identifica- 
tion of the appropriate cluster state allowing the realization 
of the DJ algorithm in the present work: Fig. [Tta) shows the 
graph corresponding to a two-dimensional six-qubit cluster 
state where the numbered vertices stand for physical qubits. 
These qubits are equally distributed among two photons, la- 
beled as A and B: qubits 1, 2 and 3 belong to photon A 
and interact by two controlled-Z gates represented by ver- 
tical connections on the graph, while qubits 4, 5 and 6 are 
associated to photon B. As usual in the one-way model, 
it can be useful to think of the distinct horizontal qubits as 
"the original [logical] qubit at different times" lfl6Tl : indeed, 
we identify the logical qubits 1 and 2 with physical qubits 
1, 4 and 3, 6, respectively. The ancillary qubit 3 is repre- 
sented by qubits 2 and 5. The "E cluster" just described is 
our quantum computer; we can show that the choice of the 
measurement sequence for the two qubits associated to the 
ancilla leads to the evaluation of both the balanced fs and 
the constant fc functions. This implies that, in the E cluster, 
qubits 2 and 5 play the role of the oracle, while the remain- 



CNOTjy indicates a controlled-NOT gate between logical qubits i and j. 
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ing qubits constitute the tools we have at our disposal to dis- 
criminate between a balanced and a constant function. To 
better understand this feature it is useful to consider the cir- 
cuit representations associated to the realization of the two 
considered global properties of the 2-bit boolean function /. 

The proposed measurement configurations are the follow- 
ing: 

1. balanced function - By measuring qubits 1, 3 and 5 in 
the bases £?i(0), -Bs(O) and -85(71") we implement, at the log- 
ical level, the two Cnot gates (Uf B ) needed to implement 
the oracle function fg (see Fig. Then we proceed with 
the measurement of the output qubits 4, 6 and 2 in the bases 
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C%' and C* 2 
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2. constant function - We measure qubits 1, 3 and 2 in the 
bases -Bi(O), #3(0) and C 2 °\ These operations implement, 
at the logical level, the identity gate Uf c needed to imple- 
ment the oracle function fc (see Fig. |3). Then we read the 
result of the computation on the output qubits 4, 6 and 5 by 
measuring them in the bases c|°\ and Bs(tt). 

We define Bj(a) = {\a+)j, \ct-)j} with \a±)j = 
^(|0) ± e-^ll),-), while cf = {|0)„ |1),} is the com- 
putational basis for the Hilbert space associated to qubit j. 
The above sequences of single-qubit measurement lead us 
to the circuits shown in Fig. |2]and[3j in particular, we can 
see the elements of the circuit realizing the unitary trans- 
formation Uf for the balanced function fs and the constant 
function fc, as well as single-qubit Pauli gates. Here and 
in the following we indicate with Z (X) the Pauli matrix a z 
(cr x ). For a given basis Bj{a), we introduce the quantity Sj 
whose value is (1) if the measurement result is equal to 
\a+)j (|a_)j) and equivalently for the cj basis. Accord- 
ing to the algorithm, we will expect, as outputs, the state 
|1 © si © s 5 )i|l © s 3 © s 5 >2|l © 55)3 for the balanced and 
I s 1 © s 2 ) 1 1 s 2 © S3) 2 1 1 © S2 ) 3 for the constant function. In the 
previous expressions we take into account the feed-forward 
corrections of the Pauli errors. 

Experimental preparation of the cluster state. - Refer- 
ring to Fig. [TJ the two-dimensional cluster state IE) is ob- 
tained from a six-qubit hyperentangled state 11711 . |HEg), 
whose graph is shown in Fig. [Tib) |l3, 171. Our experi- 
mental setup adopts a source of two-photon states based on 
a Spontaneous Parametric Down-Conversion (SPDC) pro- 
cess where the two particles are entangled at the same time 
in the polarization and in two linear momentum degrees of 
freedom (DOFs). By a proper interferometric setup JUl it 
is possible to measure the two spatial DOFs; these vari- 
ables, labeled as the "right/left" momentum (r/£) and the 
"external/internal" momentum (E/I), are both associated to 
each of the eight modes on which the two photons are emit- 
ted. A detailed description of the experimental setup, en- 
abling the transformation from the six-qubit hyperentangled 
state |HE 6 ) into a linear cluster state, can be found in re- 
cent papers jit [l3Tl . It is interesting to note that the transi- 



tion from the one-dimensional linear cluster state to the E 
cluster considered here is entirely determined by the choice 
of the controlled-er z (CZ) operations corresponding to the 
vertical links in the graph (see Fig. Ufa)). These gates are 
optically implemented between couples of qubits belonging 
to the same photon. The graph associated to the six-qubit 
cluster state |E) exhibits two links between qubits 1 and 2 
(encoded in the E/I momentum DOF and in polarization, 
respectively) and between qubits 2 and 3 (with qubit 3 en- 
coded in the r/i momentum DOF), hence the corresponding 
CZ12 and CZ23 logic gates only involve qubits belonging to 
photon A, as already noticed above. The optical implemen- 
tation of the two controlled-Z gates is realized by means of 
two half-wave plates, as shown in Fig. [TJ C )- 

In order to give an explicit expression for the six-qubit 
cluster state produced in the present experiment we point out 
that the experimental hyperentangled state, which we label 
as |HEg), does not coincide with the hyperentangled state 
|HEg) corresponding to the graph in Fig. [TJa) and instead 
satisfies the relation 

|HEe) = H4Z5H5X6H6IHE6) = 

= 4=(|00>14 + |ll)l4) © ^(100)25 - |11>26)® 



-L(|oi) 36 + |io) 3 6), 



where Hj is the Hadamard gate on qubit j and Xj (Zj) is the 
o~x (°z) gate on the corresponding qubit. For the E cluster 
represented by the graph in Fig. [TJa) we can write that 



|E) =CZ 12 CZ 23 |HE 6 ). 
Combining Eq. © with Eq. (UJ we get 

|E) = CZ 12 CZ 23 |HE 6 ) = H 4 Z 5 H 5 X 6 H 6 |E) 
1 



(2) 



■X\EE)\<t>-)M) + \EE)\<t> + )«\Zr) 
\II)\^+)^\rl) + \II)\cp-)A^)) 



(3) 



for the six-qubit two-photon E cluster state generated in 
the laboratory. In the above expression the states |(/> + ) w 
and \<fi~)iT are the two polarization Bell states. As usual 

HEED], 

we refer to |E) and |E) as the state in the "clus- 
ter" and "laboratory" basis, respectively. As shown in Fig. 
[TJc) we transform the hyperentangled state |HEe) into the 
cluster state |E) by applying two CZ operations. 

Experimental results. - In order to characterize the gen- 
erated |E) state we measured the witness operator W = 



3 - 2(IlLi ^t 1 + IlLi 22-^i±i) 13 (see Q for the 



definition of the We found (W) = -0.333 ± 0.002, 
demonstrating a genuine six-qubit entanglement |19]. Since 
it is possible to show [18, 2^1 that the fidelity F satisfies the 
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TABLE I: Experimental probabilities of the obtained output states 
for balanced and constant function, with (FF) or without (no-FF) 
feed-forward. We indicate in bold character the data corresponding 
to the expected outputs. 





/b: Balanced 


fc '■ Constant 


Output 


No-FF(%) 


FF(%) 


No-FF(%) 


FF(%) 


|000) 


0.8 ±0.1 


0.9 ±0.1 


0.7 ±0.1 


0.9 ±0.1 


|001) 


2.7 ±0.2 


2.6 ±0.1 


77.5 ± 0.5 


75.2 ± 0.2 


|010) 


1.4 ±0.2 


1.4 ±0.2 


1.2 ±0.1 


1.4 ±0.1 


ion} 


15.5 ±0.5 


14.1 ± 0.1 


3.5 ±0.2 


3.4 ±0.1 


|100) 


1.3 ±0.2 


1.0 ±0.1 


1.2 ±0.1 


1.0 ±0.1 


|101) 


2.4 ±0.2 


3.4 ±0.1 


13.5 ±0.4 


14.1 ± 0.2 


|110) 


0.4 ±0.1 


1.4 ±0.1 


0.4 ±0.1 


1.4 ±0.1 


1111) 


75.5 ±0.6 


75.2 ±0.2 


2.0 ±0.2 


2.6 ± 0.1 



relation F > |(1 — (W)), a lower bound for the fidelity is 
easily found: 

F > 0.667 ± 0.001. (4) 

Let's now turn to the DJ algorithm. We performed the 
sets of single-qubit measurements stated above and found 
the results presented in Table HI here we show the probabil- 
ities of the outputs of the computation when no Pauli errors 
are present (No-FF). This corresponds to consider only the 
case where si = S3 = S5 = for the balanced function 
and si = S2 = S3 = for the constant function. We also 
show the results obtained by considering all possible outputs 
and applying the feed-forward (FF) operations correcting the 
Pauli errors (see also Figs. |2]and[3]). It is worth noting that, 
since the output of the computation is read in the {|0), |1)} 
basis, the FF is a relabeling feed-forward, i.e. "the earlier 
measurement determines the meaning of the final readout" 
(see "Grover's search algorithm" section of ifioll or the end 
of section II in H). 

It is also important to notice that the physical qubits con- 
stituting the |E) cluster were actually measured in the appro- 
priate laboratory basis, which differs from the cluster basis 
when a single-qubit gate acts on the considered qubit; refer- 
ring to Eq. y), this corresponds to the case of qubits 4, 5 
and 6. 

The experimental results are in good agreement with the 
theoretical predictions for both functions. The main discrep- 
ancy resides on the output probabilities of the states |011) for 
Jb and |101) for fc- These states differ from the expected 
outputs in the value of the logical qubit T. This is mainly due 
to the non perfect interference visibility associated to the E/I 
momentum DOF (V ~ 70%). We attribute this to the diffi- 
culties in obtaining a perfect mode matching in the second 
interferometer (see JH for more details). 

Conclusions. - We have presented an all-optical imple- 
mentation of the DJ algorithm for 11 = 2 qubits. For this 



purpose, by taking advantage of the generation of a six-qubit 
two-photon hyperentangled state, we created a novel, high fi- 
delity, two-dimensional six-qubit cluster state that represents 
the first step for the realization of the algorithm as a one- 
way quantum computation. We were then able to evaluate 
a two-bit balanced function as well as a constant one and to 
discriminate between them in one single run of the executed 
program, in contrast to the three runs needed with a classi- 
cal computer. The correct output is identified at a frequency 
of almost 1kHz without feed-forward, a result which over- 
comes by several orders of magnitude what can be achieved 
with a six-photon cluster state, according to the current op- 
tical technology. By using all possible detection outputs and 
applying the feed-forward corrections we could obtain a fre- 
quency 8 times larger. Note that our experiment was actually 
performed with four detectors Jit]. In order to consider all 
the possible outcomes at the same time we would need 16 
detectors. 

The experimental results demonstrate the correctness of 
the proposed algorithm implementation and represent the 
first proof of such a computation with a two-bit function in 
the framework of the one-way model. 

We thank R. Jozsa for useful discussions. 
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